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We analyze the anomalous behavior in hole-doped cuprates near | doping in terms of the 
commensurate-incommensurate transition of a stripe phase. Based on an effective Ginzburg-Landau 
theory appropriate for weak pinning, we calculate the commensurate-incommensurate transition 
point and the energies of the phason and amplitudon collective modes. Using experimentally avail- 
able parameters, we estimate the phason gap (pinning frequency), the conductivity and the contri- 
bution of the phason mode to the dielectric function. 



I. INTRODUCTION 

Many of the high-T c cuprate materials show signa- 
tures of static or dynamic spin and charge modulation, or 
stripes The effect of stripes on high-temperature 

superconductivity is a subject of ongoing debate: Some 
researchers consider them central to the mechanism of su- 
perconductivity, others consider them detrimental, and 
yet others argue that stripes are irrelevant. The situ- 
ation is further complicated by the difficulty of detect- 
ing the presence of stripes, particularly when they are 
only weakly pinned to the lattice, or dynamical. The 
former case is typically associated with zero-frequency 
incommensurate-momentum response in neutron scatter- 
ing 3 , while the latter is associated with finite-frequency 
incommensurate- momentum response 0]. The difficulty 
of direct detection has motivated other approaches to 
stripe detection based on their effects, e.g., on the phonon 
modes or lattice structure [1, 0|- In particular, one 
naturally expects that in the presence of sharp stripes 
there will be a rearrangement of the collective mode spec- 
trum, with the formation of local modes split off from the 
bulk continuum, or edge modes [HQ. 

One of the most striking effects in the striped 
cuprates is the suppression of the superconductivity near 
the | hole-doping, where (based on neutron scatter- 
ing) the stripes become stationary: the so-called "|- 
anomaly" At | doping, stripes are locked in 

with the lattice resulting in a commensurate modulation 
wavenumber. This suppresses the charge and spin fluc- 
tuations, which are considered to be important for super- 
conductivity. Away from ^ doping, on the other hand, 
the favorable wavenumber for the stripe structure is ex- 
pected to be incommensurate and proportional to the 
hole concentration. Naturally, the incommensurability 
renders the stripes more dynamical. While this anomaly 
strongly suggests that the commensurate (C) stripes sup- 
press superconductivity, the question of whether incom- 
mensurate (IC) (or fluctuating) stripes promote the su- 
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perconductivity remains. We believe that study of this 
interplay is important for developing an understanding of 
the high-temperature superconductivity in cuprates. As 
a first step, in this paper, we discuss the C-IC transition 
near the commensurate doping, calculating the transi- 
tion point within a weak pinning approximation. We 
also calculate the energy of two charge-spin-lattice ex- 
citations, namely, the phason and the amplitudon. In 
the last part of the paper, we provide relations between 
Ginzburg-Landau parameters (order parameters, energy 
coefficients) and physical quantities (conductivity, dielec- 
tric function) by considering the response to an electric 
field. 



II. MODEL 

We model the stripe phase by a simultaneous unidi- 
rectional sinusoidal charge density wave (CDW) and spin 
density wave (SDW) with wavelengths A = and 2A, 
respectively. Here a is the lattice constant, G is the 
reciprocal-lattice vector, and q is the ordering vector. We 
introduce rp(xi) and x( x i) as t ne complex order param- 
eters for the CDW and the SDW, respectively. The real 
part of ip{x) gives the deviation of the charge density 
from the uniform state: 

i>(x)+^*(x)=p(x)-p , (I) 

where po is the average hole-density per site: p = z for a 
doping z. The real part of x( x ) gives the staggered spin 
density: 

X (x)+x*(x) = (-lV /a o-(x). (2) 
The periodicity of the CDW, v, is characterized by 



If the CDW and lattice are commensurate, we can write 
v = L/L with relatively prime integers L and L. We 
extend v = L/L to irrational numbers for the IC case by 
taking L — * oo and L — > oo. We define -r as a commen- 
surability; thus, the commensurability for the IC case is 



2 



zero. The order parameters have the following forms: 

ip(x,<f>) = pexp[i(qx + (4) 
X{x,4>) = -imexp[i(|a; + ^)]. (5) 

We investigate the free energy near the commensurate 
doping z c = i. Here we consider only the following two 
cases as the groundstate candidates. One is the IC case, 
where the periodicity is given by v = (L = oo) as a 
function of the doping z. The other is the simple C case 
with v = L (finite L) independent of the doping z. We 
denote q especially in the IC case by Q{z) = 2zG. 

We introduce the Ginzburg-Landau functional (per 
site): 



F 



7T fbl>(x),X(x)]- 



(6) 



In general, / can be written in the following form for a 
one-dimensional CDW with commensurability [TTL fT3| : 

f[il>(x), X (x)] = f mx),x(x)}+p(x)Re{^(x) L } 

+. 9o |[Q(z) + iV]V(x)| 2 . (7) 

Here, /o is concerned with the spin and charge ordering, 
but not with the pinning or deformation of the CDW. 
The second term in expression (0 is introduced to rep- 
resent the pinning energy of the CDW. The third term, 
which originates from the gradient term in the free en- 
ergy, gives the preference for Q{z) (g > 0). We introduce 
this term to describe the energy which arises from the 
deformation of the CDW; however, to simplify the cal- 
culation, it does not include the local effects (e.g., those 
from discommcnsurations [12j |) but only gives the prefer- 
ence among different wavevectors. We assume that /o is 
of the form 

fobP(x),x(x)] = r \x(x)\ 2 + u Q \x(x)\ 4 

-so |x(z)V(z)| +v \i>(x)\ 2 , (8) 

with the parameters chosen in such a way that the charge 
order is induced by the magnetic order [l3|. Expanding 
p(x) in harmonics, 



the pinning term disappears for incommensurate cases, 
since p L — > for L — > oo. The tq and uo terms gov- 
ern the continuous transition between the stripe and the 
uniform states at a large doping. Since we are interested 
in the stripe state, we take vq < 0, uq > 0, so > and 
vq > 0. ro(z) vanishes at the order-disorder transition 
point zq. The region between the C-IC transition points 
z c _ < z < z c+ around the C point is the C region; 
outside this region the system is in the IC phase (Fig.[T]). 
We will find in the following that the width of the C 
region rapidly shrinks with increasing L. 
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p(x) = ^2picos(IGx), 



(9) 



FIG. 1: The C-IC transition near the commensurate doping 
z - i 



Although the term p L can give large negative energy 
for large p, this does not cause any problem so long as we 
restrict pc to not too large values, so that the p^. term 
is not dominant. In other words, this term makes only 
a minor modification of the free energy, which selects 
between the C and IC equilibrium states. 

We make the following assumptions in this paper: 

1. \pi\p L ~ 2 -C vq (weak pinning). 

2. Sq <C uqVq (weak coupling). 

3. ro is constant around the C region (z c ± -C zq). 
I m< 0.1. 

5. A(= -^j) < 1 (equivalent to s J$ ^o)- 

Assumption 4 is consistent with the neutron scattering 
experiments where it is roughly estimated that 2m ps 0.1 
for striped cuprates [14j . In assumption 5, the inequality 
sign stands for "to be on the order of or less than" , and 
the assumption reflects the fact that the charge order is 
much weaker than the magnetic one in metallic stripes 
due to the strong fluctuations. 



the free energy per site can be rewritten as 

F{p,m,<l>) = F ( -°'>(p,m)+p 1 p L 
+9o[Q(z) - q} 2 p 2 , 
F^ Q \p 1 m) = tqui 2 + ugm 4 — s$m 2 p + vop 2 

where 

Pi = T^Pi cos(0 o ), 9o = ^9o- 



(10) 
(11) 

(12) 



For the equilibrium state, 4>q is taken so that p\ cos(0o) = 
— \p\\ (i.e., pi is always negative, and the sign of p± de- 
termines the location of the stripe center). Notice that 



III. RELATIVE STABILITY OF C AND IC 
PHASES 



Now we evaluate the order parameters for the C and IC 
cases so that p and m minimize the free energy. First, in 
the IC case, it follows from |^ = 9£ = and Fin = F^ 

1 am op 

that 
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(13) 
(14) 
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The following condition is required for m 2 c > and 
Pic > 0: 



Sq < 4u v , 



(15) 



which is satisfied by assumption 2. Substituting Eqs. (|13p 
and (|14p into Eq. (| 1 1 1) , we obtain the condensation energy 



(16) 



Next, in the C case, the magnetic order parameter is 

(17) 



2 M + S pc 2 

m n = m m 



2u 



ic- 



The charge order parameter for the C case satisfies the 
equation 



Pc 



L 1 '-Ap c +B = : 



(18) 



where 



2v ~^ + S~g G 2 iz-S,) 
A = ^2 — ; i V (19) 



L\Pi 



D 



Vo\ so 
2m L \pi 



(20) 



Since £? > and A < 0, the following condition is re- 
quired for the existence of a positive root: 



B-(L-2) 



\A\ 
L-l 



L-l 
L-2 



< 0. 



(21) 



Under this condition, Eq. (fTS)) has two positive roots. 
The smaller one corresponds to the minimal energy in 
the C case, while the other corresponds to the maximal 
one; therefore pc is given by the smaller positive root of 
Eq. (18!. 

A is a function of z and L. For later use, we define Aq 
and Bq independent of z and L near the commensurate 
point: 



A = L\pi\ A + 2Z 2 = 
ko|so 



AuqVq - g| 

2u 



B = L\p 1 \B = 



2un 



where 



4g G 2 (z-^L\ 



2GJ 



(22) 
(23) 

(24) 



Note that, near Bq depends on z through ro, but 
not so far below zo (from assumption 3). From now on, 
we also use the notation Z with a subscript to represent 
that corresponding to z with the same subscript through 
Eq. (24]). 

We discuss the C-IC transition as a function of the 
hole concentration z. We will evaluate the difference of 
the free energy between C and IC below. However, F\q 



is not well-defined at Z = since there is no IC state 
at Z = 0. To avoid problems with this discontinuity, we 
first define 



Pic 



Bp 
\A \ 



\r \s 



4u v 
Fic = F(°\p IC ) 



(25) 
(26) 



as the order parameter and the IC free energy for any Z 
including Z = 0, and then we take the commensurability 
at the Z = point into account after comparing the 
continuous free-energy functions if needed. 

The difference between the C and IC free energies is 
[from Eqs. (fl3j). (|TT|). and (fT8])] 



AF = Fic - F c = AF (0) 

2 \ L 
1 
L 



+ 1 



Ipxl^-Z 2 ^ (27) 
1 + Z 2 ) p 2 c 



\Ap\ picpc 



-Pio 



(28) 



where Z 2 = ^ 



The roots of AF = for L > 3 are 



(Ai^O) _ 

Pc - 



L-l±^l-L(L-2)Z 2 
(L-2)(l + Z 2 ) 



-Pic- 



(29) 



Here the smaller root gives the C-IC transition, because 
should be positive so that d ^f^) < at the tran- 



sition point. We require Z 2 < 



l 

L(L-2) 



so that the C-IC 



transition exists. For L > 3, p c F has only one mini- 
mum at Z = 0. We find that pc is larger than p\Q at Z c 
(< Z ): 



(AF=Q) 

Pc 



PIC 



1 



L-2 ~. 



> Pic- 



From Eq. (|18[) . we can write 
Pic 



Pc 



L-l 

Pc 



(l + Z 2 ) \A\ 



(30) 



(31) 



pc and pic vanish together at the order-disorder transi- 
tion point Z = Zq. Note that we take pic as a constant 
for Z far from Zo, but it is a function of Z near Zq which 
vanishes at Z = Zq. The order parameter at the C- 
IC phase transition, pic(Z c ) 0), should satisfy both 
Eqs. (|30 | and (|3T |) . 

For Z* < 1, from Eq. (pj) . 



JAF=0) 

Pc 



L-l 



L-l 

Pic 



1 



(L-2)(L-l) g2 



(32) 



Substituting this and Eq. ([3D]) into Eq. ([3T|) and keeping 
terms up to the order in Z 2 , 



1 



L-2 



-zi 



(1 + LJT) 



L 2 (L-3) 



X - 1 



2? 



4 



2X 



1- L(L-3)X' 



(33) 



where X = fc*M 



Since X 2 ~ from assumptions 1 and 2, Z 2 (= 
is represented by a simple function of \pi \ : 



Z. 



2pL -2 \Pl\ 



(34) 



By virtue of assumption 1, it is apparent that Z^ <C 
1, which is consistent with the case we consider here; 
therefore it is confirmed that Z c in Eq. (|34|) is the C-IC 
transition point for the case where assumptions 1 and 2 
hold. 

The width of the commensurate region is given by 




(35) 



The commensurate region becomes exponentially narrow 
as the CDW period L increases. Therefore, only com- 
mensurate stripes with the relatively small L may exist 
in practice. Note that, here, we do not consider the case 
where discommcnsurations exist II 211 . 



At a commensurate point, the condensation energy is 
given by 



E con - -E$ n - Ec, 
41 = l\ro\m 2 c , 

Ec = ^-\pi\ P L - 



(36) 
(37) 

(38) 



IV. COLLECTIVE MODES: PHASONS AND 
AMPLITUDONS 

Now we consider the slightly displaced CDW to study 
collective excitations around the equilibrium state. First, 
we discuss the phason and then the amplitudon. 

For the phason, we introduce an ^-dependence into the 
phase of the order parameters: 



pexp 



i(qx 



x{x,(j)) — — lmexp 
tp(x) = £cos(fcx + cf>i), 



-;■<) 
L 



<£<£)■ 

L 



H 2 X ' 2L 



(39) 

(40) 
(41) 



where k = with N > 2. The free energy is given by 



Fo(p) 



tlx 



y < Mi>(x), X (x)] +p(x)Re{^(x) L } +g 



Q(z) 



kt; sm(kx - 



W(x)\< 



F (p) 



F^(p)-\p 1 \p L +g [Q(z)-q] 2 p 2 , 



1 ,- , L 1 

- 2 \pAp l - 



5 k ,o cos(20i) ( 1 g p 2 fc2 _ 1 + <5 M cos(20i) p 2 



2 L 2 



L 2 



L 2 

y\Pi\p 



L-2 



gok 2 



(42) 

(43) 
(44) 

(45) 



Assuming that p depends on time t only through £(t), 
the kinetic energy is given by 



T 



f dx 




1 v 


dt 



h P 2 l + 5 kfi cos{2(j) 1 ) fd£ 
L 2 2 I dt 



(46) 
(47) 



where to is a Ginzburg-Landau parameter of the kinetic 
energy and defined by this equation. Thus, the La- 
grangian corresponding to the fluctuating CDW is 



C = \ [l + 4,o cos(20i)] 



top" 



d 4) 2 -h{^\p L - 2 + ~^ 2 )e 



const. 



(48) 



Solving the Euler-Lagrange equation of motion, 



oc cos(fi<^ + <j> 2 ), 



i 

to" 



■\vi\p 



L-2 



9ak 



(49) 
(50) 



f2 v is gapless in the IC case and gapped in the C case. 



From Eq. (j3Tj) . the gap at k = can be evaluated as 
L 2 



A 



L-2 



(51) 



— \PMc 2 [1 + (L ~ 2)X] [l — (L — 2)Z 2 } , 



2t 



(52) 



where we have used Z 4 < Z^ <C 1 in the commensurate 
region [see Eq. (|34p]. In Fig. [2J the doping dependence 
of the phason gap is shown for the L = 4 case. 
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FIG. 2: The phason gap around g doping (L = 4). The width 
of the commensurate region is z c + — z c 



In Appendix [B] to is roughly estimated as ~ p c 2 x 
10~ 32 [eV s 2 ] using a single-band Peierls-Hubbard model. 



where Eq (— |pi|p 4 ) is the commensurability energy. Eq 
is roughly estimated as U 7 m 8 x 10~ 2 [eV] in Appendix iBl 
(U is the on-site Coulomb interaction in the single-band 
Peierls-Hubbard model). Therefore, the order of the pha- 
son gap at z — i is 



Ms) 



(7 2 777 X 



10 15 [Hz 



'U^m 4 [eVl 



(54) 



We expect that U < 10 [eV] and m < 0.1 (assumption 4), 
and thus find the order of A v (|) is less than 100 [meV]. 
For a smaller U , it can be even much smaller than 100 
[meV]: e.g. for U » 1 [eV], we find A v (±) < 1 [meV]. 

Next we consider the amplitudon mode. For this pur- 
pose, we write the order parameters in the following form: 

i/>(x, 4>) = p [1 + A(x)} exp + %)] , (55) 

X (x, <f) = -lm [1 + ^(x)] exp U\x + (56) 
^(x) = Ccos(fca; + ^i). (57) 



x 10 16 [Hz], 



(53) Then, the free energy is given by 



F(p, = f y {foMx), x(x)} + p(x)Re {^{x) L } + g [Q(z) - qf \^(x)\ 2 + g Q H sm(kx)} 2 p 2 ) 



F (p)+F^\p)C 2 , 
1 + 4,ocos(2</>i) J |r | 2 |r |s 



U 



+ p+fvo+golQW-qfjp 2 + 



1 1 l(jL ^IpiIp l + ^V 



(58) 
(59) 
(60) 



Thus, we find the frequency of the amplitudon fluctuation 

as 



1 



M*0 = t 

to 



For -2 
— p 

UQ 



Fo so -i 

P 

it 



+ vo+5o[Q0?) 

+ ^ L \Pi\p L 2 +9ok 2 

A 2 A + ^k 2 
to 



(61) 
(62) 



The amplitudon has a gap even in the IC case (unlike the 
phason). The amplitudon gap at h doping is 



A. 



'2|ro|m 2 + nop 2 



top 2 



~ V4°l x 10 16 [Hz] 
< Urn x 10 16 [Hz], 



(63) 

(64) 
(65) 



where U is in the unit of eV, and U <C 10 [eV] (see 
Appendix [B| . 



So far, we have discussed the one-dimensional collec- 
tive modes, that is the modes that are uniform in the 
direction along stripes. Next we consider the meander- 
ing mode of stripes, i.e. the transverse mode with the 
wavevector oriented along the stripe direction. This case 
is a simple extension of that for the one-dimensional pha- 
son. The only difference is that now we impose a y- 
dependence on £ and introduce the energy cost associ- 
ated with the spatial variation in the y-direction. Then, 
we obtain the following Lagrangian 



C = 



i[l + ^ cos(2^)]^- 



dy_ 

V v 



d£\ 2 _ho (dC 2 
dt ) t \dy 

L-2 



+ 9ok 2 x £ 



const. 



(66) 



where /iq is the energy coefficient of the ^-derivative term. 
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Solving the Euler-Lagrange equation of motion, we find 

(67) 



£{y,t) oc cos(k y y + fl m t), 
2 1 



to 



g k x 



(68) 



In a similar manner, we find the two-dimensional am- 
plitudon as 



C(y,t) oc cos(k y y + Q A t), 



hok 2 + g k 2 



A3, 



(69) 
(70) 



V. RESPONSE TO ELECTRIC FIELD: 
CONDUCTIVITY AND DIELECTRIC 
FUNCTION 

To investigate collective oscillation of stripes which can 
be experimentally examined, we consider the meandering 
mode contribution to the conductivity. Since we do not 
include impurity pinning effects in this paper, our treat- 
ment is equivalent to the classical particle model in the 
one-dimensional case. In the presence of an electric field 
E = Eg ex.p[i(ky+u>t)] along the x direction, the equation 
of motion is given by (see Appendix [Cl or reference (l5j ) 



d 2 Z(y,t) , d£{y,t) h d 2 £,{y,t) 



dt 2 



dt 



t dy 2 



= — E Q exp[i(fcy + u^Tl) 



where /i is the effective mass of the CDW (see Ap- 
pendix [B]). The particular solution of this equation is 



£(!/,*) = 



ep G 



E 



n [n m (o,k)} 2 - lu 2 +iojV 

Therefore, an induced current-density is 



(72) 



. e(p /a 2 c) dj 1 iuu 2 E(y,t) 

Jm{y,t) = -p, — = —7^2 2 1 .. ,-n ' ( 73 ) 



G dt 4tt fll - lu 2 + iljT ■ 



where 



2 = lW_ m _ 3 , 

a cfi 



(74) 



and c is the lattice constant in the z-direction. 

At I doping, from the effective mass (fi « 9.2 x 10~ 31 
[kg]) estimated in Appendix [Bl it follows 




10 15 [Hz] ~ 10 14 [Hz] 



(75) 



Once we know io v accurately, the charge order parameter 
can be determined. The value of lu p will be determined 
from the phason contribution to the dielectric function 
and the phason gap, A v (see below). 



The conductivity is given by 



u(/c, to) 



dE 4tt fll - lu 2 + iuT 



(76) 



The contribution of the meandering mode to the dielec- 
tric function is 



e m (k,w) = —a{k,uj) = — v — -. 

lui lli, — lu z + livL 



(77) 



The contribution of the phason mode to the static dielec- 
tric constant is 



w£ _ 367T/9§e 2 x 10 9 
2" 



c >-o)-^_ {aq?cL2Ec 



(78) 



where Ec is in unit of hertz. From the experimental 
data for e v (u} — 0), Eq can be determined. We can 
also determine ui p from e v (u> = 0) and A^,(|), and the 
effective mass can be estimated from lu p (this would be 
close to that of a free electron). 

At g doping, it follows from the estimated phason gap 
(l54)l and plasma frequency of the CDW ([75]) that 



(79) 



eJu = 0) ~ U~ 7 x m - & x 10 _1 , 



where U is in unit of eV. 



VI. SUMMARY 

In this paper, we investigated the competition between 
C and IC stripe phases by calculating the C-IC transition 
point in the weak-pinning approximation. This allowed 
us to estimate the strength of the lock-in effect, and eval- 
uate the phason and amplitudon fluctuation frequencies, 
the ac conductivity, and the dielectric function. 

Our calculation suggests a narrow C region around | 
doping. In this region, there is a weakly-pinned phason 
mode (with a charge-order-dependent gap at k = 0). The 
phason gap at k — for g doping is estimated to be of 
the order of 10 meV. The pinning frequency increases 
with |fc|, quadratically near |fc| = and linearly for large 
|fc|. This weakly-pinned behavior should be observable 
by electronic transport experiments. In the IC regime 
(away from | doping), the phason frequency is linear in 
|fc| with no gap at k — 0. The amplitudon shows a similar 
^-dependent frequency as the phason. The amplitudon 
gap at k — in the C phase is larger than that of the 
phason, and remains finite in the IC case. 

In the presence of an electric field across the stripes, 
phason or meandering fluctuations can be induced. From 
measurements of the static dielectric function or conduc- 
tivity, it should be possible to extract the amplitude of 
the charge order parameter. These and other measure- 
ments would help constrain parameters in the Ginzburg- 
Landau model considered here. Also, direct real-space 
calculations [161 ] can decrease the number of independent 
parameters. 
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A similar model can be applied to describe 
commensurate-incommensurate transitions and response 
functions in other related classes of materials, including 
nickelates and manganites. 

This work was supported by the U.S. DOE. 



APPENDIX A: SINGLE-BAND 
PEIERLS-HUBBARD HAMILTONIAN 

Here, we consider a single-band Peierls-Hubbard 
Hamiltonian and estimate the electron-lattice coupling 
strengths. The Hamiltonian reads 

n = - ^ [ty + 47^- {wi + Wj )J c\ a c ja 

+ U E n «T n a + 4 7 E W i C io C ia 
i i,cr 



where 



Wi = ^ («<«,,»» - Ui x -i,i v +v ix , iy - v ix! i y -i), (A2) 

(c) is the hole creation (annihilation) operator, is 
the Cu-Cu hopping between the nearest neighbors (t) or 
the next-nearest neighbors (£'), and Ui and v\ are the 
displacements of the x-Oxygen and j/-Oxygen at site i, 
respectively. 

To estimate the electron-lattice coupling strength in 
the single-band model, we start from the three-band (Cu- 
O) model. The Cu on-site energy and the Cu-Cu hopping 
integral to the zeroth order in t p d are given by 

ef > = -A pd = -A + 4/3u>i, e\f = 0. (A3) 

The second-order perturbation to the Cu on-site energy 
and the Cu-Cu hopping integral are given by 



,(2) 



E 



A pd (i) 



(A4) 



A=l,2,3,4 

c (2) 1 (tpd(i,\j)tpd(j,^ji) , t pc [(j,Aji)tpd(i,\ij) 



Apd(j) 



(A5) 



where the Cu-0 hopping and d-p energy difference de- 
pend linearly on lattice displacements, 

t pd (i,\) =t pd =fau(i,X), A pd = A - 4/3w l; (A6) 

and A = 1, 2, 3, 4 refer to the oxygen ion to the right of, 
above, to the left of, and below the Cu ion, respectively; 
and Xij to that between the i and j Cu ions. By re- 
taining up to linear terms of u and v, e- and are 



approximately 



(2) ^ u ld ~ Stpdawt 
€l ~ A - 4/3^ 

Kd ( 16 ^ 8at pd \ 
A \ A 2 A J 



(A7) 



1 



+ 



t 2 

L pd 



2 \ A - 4/3wi A - 4:/3w 



t%d . 2pt 



A A 2 



(A8) 



Combining now the zeroth and second-order results, 
for the effective one-band model, we obtain 



A 4 & 4 

" Ao " a7 4 

47 ^ + const. 

mid 



4/3- 



16f3t 2 



pel 



A 2 



8at pd \ 
A J 



(A9) 



pa ' pa t . \ 

A A^-( W *+ W j) 



/ 

A 

i-4 7l ',> 



A 2 
^0 



(A10) 



with the reduced electron-lattice coupling strengths 7 
and 7': 



2at 



7 = 



pd 



mid 
a 2 , : 



7i 



/3i 2 



pd! 



2A 2 ' 



(All) 
(A12) 



Using the three-band parameters t pd = 1.3 [eV], A = 
3.6 [eV], < a < 4.5jeV/A], < (3 < 1 [eV/A] and 
K = 40 [eV/A 2 ] [H[l3j, we obtain the single-band pa- 
rameters: t = 0.47 [eV], < 7 < 3.73 [eV/A] and < 
7' < 0.07 [eV/A]. Here the value of the Cu-Cu hopping 
energy is similar to that estimated in Ref. [l8| (t = 0.43 
[eV]); therefore, it is also expected that t' « —0.07 [eV] 
from Ref. [H. 



APPENDIX B: ESTIMATION OF THE 
GINZBURG-LANDAU PARAMETERS 

Here we represent the Ginzburg-Landau parameters 
for kinetic and pinning energies, t and pi, by using the 
order parameters and other parameters of a single-band 
Peierls-Hubbard model, and estimate the commensura- 
bility energy and the effective mass of the CDW. 

We start with the single-band Hamiltonian in Ap- 
pendix [A] We can adopt a particular set of parame- 
ters obtained there: e.g., t = 0.47 [eV], t' = -0.07 [eV], 
7 = 2.0 [eV/A], 7' = [eV/A] and K = 40 [eV/A 2 ], etc. 
However, we will not necessarily require all of these val- 
ues below, because it is not our purpose to calculate the 
order parameters here. Moreover, the order parameters 
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are assumed to be p < 10~ 2 and m < 10 _1 in this paper 
{assumptions 4 and 5). Our purpose here is to estimate 
the commensurability energy, the kinetic energy, and the 
effective mass of the CDW. In the estimation of the com- 
mensurability energy, we will find that the commensura- 
bility energy directly depends on only t, U and m as long 
as the order parameters are small (p -C m, < t/U). In the 
estimation of the kinetic energy and the effective mass, 
we will find that those values depend on 7, K and p as 
long as (pj^) 2 is small (< 10~ 5 [A]). Therefore, we take 
the parameters as t — 0.47 eV, U <C 20t, and < 0.1 
[A] ; and we assume that the calculation with the result- 
ing parameters yields order parameters which satisfy as- 
sumptions 4 and 5. We also consider the terms of 7' and 
t' to be small enough to be neglected, compared to the t 
term. 

Defining the fluctuations: 



n = n—{n), u = u — \u), 
V = v — (v) , w = w — (w) , 



(Bl) 



we rewrite the Hamiltonian: 

H = Hmf + H z + 'Hint + T~i P h + H c , (B2) 

where 



Hmf = ^2 c ta c j<* ~ *y + [ u + 4 7<Wi) J 6 i;j 



(B3) 



Ui)Ui + (Vi)Vi 



'Hint — 

i 



47 ^ 

(7 

* , „,2 



(B4) 
(B5) 

(B6) 



w,;) 2 + (vi) 2 



(B7) 



According to the mean-field approximation, we neglect 
the TLi n t and H p h terms. To consider the equilibrium 
states, the terms linear in Ui and w, are set to zero, TL Z = 
0. 

We assume the order parameters to be described as in 
Eqs. (gjl and (0): 

(n ir ) + {riii) = po + bPixi, <!>) +tP*( Xi , (/>)], (B8) 
(rki) - (nil) = \x( x ii<t>) + X*( x iA)] (-l) - "^ 



Hence, (rijo-) can be expanded as 



xfa.^e 1 ** (B9) 



f + (n QCT )e'(«+*)- + (n_ QCT )e-(Q+*> 



<"-2Q CT ) = 

where Q = i. 



+ (niQvje 
( n -Ql)* 
(n-2Q<r) 



\2Qxi 



-ilQxi 



n-2Qa)e 

1 ±_ 

( n -Ql) = 7T mcl ' 2L 



2i 



1 



From H z = 0, 



(BIO) 
(Bll) 

(B12) 



(B13) 



= \^(xi+U 



i^r^Oi, <j))e lQa sin(Qa) + c.c 
K 



^(^,0)]+c.c. (B14) 
(B15) 



We substitute the Q-represented form of (n ia ) and (v,i) 
into Hmf, and expand c and as 



(B16) 



where I runs from to 2L — 1 for a commensurability , 
the sum with respect to ko is taken within the reduced 
Brillouin zone spanned by Q, and Nk denotes the number 
of k points in the original Brillouin zone. Hence, Hmf 
can be rewritten as 



H 



MF — 



J2 [< k o + iQ)<i a+lQa c ko+lQa 



+ (l+2)Q a C k +lQ a + n - C 



where 



e(k) — —2t[cos(k x a) + cos(kya)] 
—At' cos(k x a) cos(fcj,a) + U 



Pa 



U = 



U-%- sin 2 (Qa) 
K 



(B17) 

(B18) 
(B19) 



Only the states near (±f , ±5) at the Fermi level con- 
tribute to the condensation energy The energy gap 
around (±^,±-^) is approximately given by 



2A 2 



2A 2L 

JJ2L-2 



cos(<f>), 



(B20) 



where A 



2^ and D ~ 2i. The approximation here 



is justified when D ^> A and mf/ ^> /?[/. This is the 
case for L = 4 and {/ -C 20t under assumptions 4 and 5. 
Therefore, the commensurability effect on the condensa- 
tion energy, 5c, is given by 



51 



2A 



2L 



£)2L-2 ' 



(B21) 
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Now we consider the energy to be calculated from the 
meanficld Hamiltonian 



H — H-MF + He 



(B22) 



In the energy, the component proportional to m 2 con- 
tributes to the formation of the ordered state as well as 
the ro term in the Ginzburg-Landau free energy 
The condensation energy per site is given by 



E C on(4>) 



2A 2 

IT 



(B23) 



where the first term is from 7i c , and the second from 
7~Lmf', and 5 depends on the model system parameters 
(bandwidth, interactions, etc.) and doping. Note that 
Econi^) should have two components which are, respec- 
tively, proportional to m 2 and m 2L . The former is the 
^-independent condensation energy, —E^>n, and the lat- 
ter is the commensurability energy, — Eg- Setting <p = 

Ecoli can be written by 



E£l = 2A 2 



1 1 

3 ~ U 



\r \m 



The commensurability energy is given by 

Ec — E con (~) — E con (0) = 5 1 5 2 
'U' 



(U + \r \) 



D 



2L-2 2L 

in 



2 2L-1 



(B24) 

(B25) 
(B26) 



For an insulating CDW in the nearly-half-filled one- 
dimensional system, the condensation energy is roughly 
estimated as ~ E 2 /t, so that we expect 



l'-o 



U 2 

<^~U 2 [eV] 



(B27) 



for the metallic CDW in a two-dimensional system. Since 
we consider U <C 10 [eV], U + should be on the order 
of U [eV]. In this case, the commensurability energy per 
site is obtained by 



Ec 



U 7 m s x 1CT 2 [eV]. 



(B28) 



Then, we obtain the Ginzburg-Landau parameter for pin- 
ning energy 



\pi\ ~ U 7 A~ 4 x 1CT 2 [eV], 



(B29) 



where A = -£r. 

Next, to estimate the kinetic energy, we consider the 
time-dependence of the order parameter through 



(f>(x,t) = 4> + £(t) cos(kx + cj>x). 
The displacements arc written as 



Ui (t) = Ui{t)+U*{t), 

Ui{t) = 1 ^V(^^)e 19o/2 sin(^ 



dt 



J)^(f) 



dip 



k) sin2 (y) 9 2 p 2 «dCOs 2 (fca, 



(B30) 

(B31) 
(B32) 

(B33) 
l),(B34) 



where v d = ^j^- The mean square velocity of the ion 
and the total kinetic energy are 



dUi 



dt 



N ^ 



dU t 



dt 



27 \~ ■ 2 fia\ 

k sm It ) q 



222 

P v d 



1 + cos(20i)(5fc,i 



(B35) 



1 2 1 2 1„ 
^Po^v d = -p m e v d + -M 



dUi 



dt 



= ^Pom e v d 



M f2-f\ 2 . 2 (qa\ 2 2 1 + cos(2<£i)<5 fc 



(B36) 



where M, m e , and p are the ionic mass per site, the 
electron mass, and the effective mass of a hole in the 
CDW, respectively. 

On the other hand, the kinetic energy in the Ginzburg- 
Landau expansion is written by 



1 + 4,ocos(2</>i) 



2 2 1 + 4,ocos(2^i) 
to + nq p 



(B38) 
v% (B39) 



f dx 


dip 


1 ~v 


dt 



(B37) 



where both tq and t\ are independent of p. Near \ dop- 
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ing, To and T\ are estimated as 
1 



field (discussed in Sec.fV}. We consider a plane electro- 
magnetic wave with the electric field oriented in the x-y 



. « 2.9 P0 x 10- 32 [cV s 2 /A 2 ] (B40) P^e, E = E (k y /k, -k x /k) e '(^+^)— *. Thc 



M /2 7 V . 2 fqa 
Ti = — — sin — 
1 2 \K ) \ 2 



1.7 ( ^) x HT 27 [eV s 2 ] 



(B41) 
(B42) 



Therefore, the Ginzburg-Landau parameter for the ki- 
netic energy, to, is given by (for k = and (f>\ = 0) 



33.7p + 1.7 p 2 (#) x 10 5 



4.2p" 2 x 10 



-32 



x 10~ 32 [eV s 2 ] 

(B43) 



where we use -3^ < 10 1 and the assumptions 4 & n d <5- 
The effective mass is 

/i sa 2T0/P0 + ngV [1 + cos(20i)] /po (B44) 

In the case of fe = and <pi = 0, 

"9.2 Po + 4V(^) 2 x 10 4 



Po 



x 10~ 3i [kg] 



9.2 x 10" 31 [kg] 



(B45) 



For p < 10 , the electronic energy is dominant in the 
total kinetic energy. 



APPENDIX C: EQUATION OF MOTION FOR 
CDW IN ELECTRIC FIELD 

Here, we derive the equation of motion for a meander- 
ing CDW oscillating in the x-direction in an ac electric 



sponding vector potential is A = E/(— uo). 

The vector potential couples to the electrical current 
j, which can be expressed in terms of the carrier density 
and the velocity. The velocity of the density wave can 
be expressed via the time derivative of this phase, v x = 
(1/G)(p. Thus the form of the phase that couples to the 
plane-wave electric field as above is 



¥>{x,y,t) = £(y,t)e 



(CI) 



Here we take k x — > 0. Thus, the additional contribu- 
tion to the Lagrangian due to the interaction with the 
external electromagnetic field is A • j — ^i" G A x ^. Com- 
bining it with the bare Lagrangian, Eq. (|66p . we obtain 



C 



t p 2 f dy 



<i£ \ 2 ho ( <i£ x 2 



dt 



t \dy 



dy epo E x 



V y G 



-lLj) 



(C2) 



Adding a dissipation term, we obtain the equation of 
motion: 



t dy* 



dt 2 



dt 



eG 



E x (y,t). (C3) 
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